Dynamical control of quantum state transfer within hybrid open systems 
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We analyze quantum state-transfer optimization within hybrid open systems, from a "noisy" 
(write-in) qubit to its "quiet" counterpart (storage qubit). Intriguing interplay is revealed between 
our ability to avoid bath-induced errors that profoundly depend on the bath-memory time and 
the limitations imposed by leakage out of the operational subspace. Counterintuitively, under no 
circumstances is the fastest transfer optimal (for a given transfer energy). 
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Commonly, manipulations of quantum information can 
be schematically divided into three stages: "writing-in" , 
"storage", and "reading-out" [![. Realistically, some 
systems are better suited for writing-in or reading-out 
than for storage, and vice versa. This has prompted 
the suggestion of hybrid, composite quantum systems @- 
[Bj]: quantum operations are rapidly performed and effi- 
ciently written in a qubit susceptible to decoherence, e.g., 
a Josephson (superconducting) qubit; then the quantum 
information is transferred (directly or via a link) to a 
storage qubit resilient to decoherence (encoded in an en- 
semble of, e.g., ultracold atoms); then, on demand, trans- 
ferred back and read-out from the same fragile qubit. 
Our aim is to examine a strategy for maximizing the av- 
erage fidelity of quantum state-transfer in such hybrids, 
from a subspace fragile against decoherence to a robust 
subspace, by choosing an appropriate dynamical control 
field. 

To this end, we resort to a novel general approach to 
the control of arbitrary multidimensional quantum op- 
erations in open systems described by the reduced den- 
sity matrix p(t): if the desired operation is disturbed 
by linear couplings to a bath, via operators S B 
(where S is the traceless system operator, and B is the 
bath operator), one can choose controls to maximize 
the operation fidelity according to the following recipe, 
which holds to second order in the system-bath cou- 
pling (Suppl. Info.): (i) The control (modulation) trans- 
forms the system-bath coupling operators to the time- 
dependent form S(t) £5 B(t) in the interaction picture, 
via the rotation matrix Ei(t): a set of time dependent co- 
efficients in the operator basis at (Pauli matrices in the 
case of a qubit), such that: 

S(f)=5>(i)^-. (1) 
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(ii) This allows to write the time-independent score ma- 
trix, describing how the fidelity scores (changes) for each 
pair of basis operators applied: 



r y = M[tr i ,a^)M]M'>, (2) 



where the overline is an average over all possible initial 
states, (iii) Using By one arrives at a simple expression 
for the average fidelity of any desired operation (within 
the stipulated second-order accuracy): 

/oo 
dwG(ui)F(t,uj), 
-OO (3) 

F(t,u) =t- 1 e t , < (w)r y e t * J .(w), 

where e t ^(uf) is the finite-time Fourier-transform of the 
rotation matrix £i(t), and the coupling spectrum G(us) 
is the Fourier-transform of the bath memory (correla- 
tion) function (j3e l ^ Bt Be^ l ^ Bt ^. Namely, the modula- 
tion (control) spectrum F(t, uj) is defined according to the 
operation, via the ITy score matrix, (iv) This fidelity is 
maximized by the variational Euler-Lagrange method [6] , 
which minimizes the overlap between G(u>) and F(t, u>) 
under the constraint of a given control energy or action. 

We use this general approach to optimize a reliable 
transfer of a quantum state from a fragile qubit to a 
robust qubit. We choose to focus on the case of two reso- 
nant qubits with temporally controlled coupling strength. 
The free Hamiltonian without decoherence is then 

where H c (t) is the Hamiltonian for the controlled inter- 
action between the qubits, V(t) describing the ajustablc 
amplitude of the interaction (see Fig. [1} inset, for an ex- 
ample where the interaction amplitude is ajustable using 
an external laser field, as in |2j). The system-bath inter- 
action Hamiltonian is taken to be 

Hi = S®B(t)=&W ®B(t), (5) 

where B(t) is the bath operator B rotating with the free 
bath Hamiltonian Hb ■ This model represents proper de- 
phasing in the source qubit 1 due to the bath operator B, 
whereas the target qubit 2 is robust against decoherence. 
This model can be generalized to any degree of asymme- 
try between the decoherence properties of the two qubits. 
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Equations (|4j)- ((5j) conserve the parity of the number 
of excitations. Hence, the full two-qubit system can be 
split into two subsystems O = spandg^} , ei<7 2 )} an d 
£ = spanjl^!^) , le^)}, O and E standing for odd and 
even excitation numbers, respectively: 



Hi —Hq + H& 
H =V(t)&? 



(6) 



£ 



B(t) 



where the appropriate Pauli matrices in the O {£) 
subsystems are: a® = \gie%) (e.\gi\ + H.C., af = 
\e1g2) (ei5 2 | - \91e2) (gie2\, of = 1.91.92} (eie 2 | +H.C. and 
a% = \e\e2) (eie2| — 1 51(72) (<7i<72|- In essence we have one 
resonant and one non-resonant two-level system, both 
coupled to the same dephasing bath, which renders them 
inseparable. Both are subject to the same a x control, 
which must be chosen to maximize the fidelity of a rota- 
tion in the O subsystem while keeping the £ subsystem 
unchanged. 

The accumulated phase 



V{t')dt' 



(7) 



is our control function. In the ideal case, without de- 
coherence or leakage, the state transfer from qubit 1 to 
qubit 2 can be perfectly realized if at the final time, tf, 
the phase <p(i) satisfies 4>(tf) — \, whence any initial 
state of qubit 1 is mapped onto that of qubit 2 (initially 
in the ground state) 

(a\gi) + /3|ei» \g 2 ) -> \gi) (a\g 2 ) - i/3|e 2 » , (8) 

for any normalized a, j3. Here the states \g\) (|g 2 )) and 
ei) (l e 2)) are respectively the ground and the excited 
states of the source qubit 1 and the target qubit 2. 

There are two conflicting noise (error) considerations 
for the transfer, each affecting a different subsystem: (i) 
In the presence of interaction between the source qubit 1 
and the bath, the longer the information stays in qubit 1 
the lower the fidelity of the transfer (manifest in subsys- 
tem O). (ii) On the other hand, if we make the transfer 
extremely fast, it may result in population from \g{) \g 2 ) 
leaking into \e\) lea), thus lowering the fidelity of transfer 
(manifested in subsystem £). Such leakage [(| 0] signifies 
the violation of the rotating wave approximation (RWA). 
Namely, fast modulation V(t) may incur unwanted, off- 
resonant transitions if the transfer rate is comparable to 
the energy difference (level distance) of the qubits, ujq. 

We first focus on bath-related errors (i), assuming that 
the RWA is valid, i.e., there is no leakage because of the 
RWA violation. This may be the case if the transfer time 
is much slower than the energy separation ujq. This is 
also true when the non-RWA terms simply do not exist, 
such as in 2D or 3D Heisenberg interactions (of the form 



® ai 2) + a { y ] ® a y 2) or ai 1] ® + ar y 1] <g> a y 2) + &P ® 
a z , respectively) where only number-conserving terms 
exist. The control Hamiltonian Hc(t) then has the RWA 
form 0: 

H c (t) =V(t) fleisa) (<?ie 2 | + | 9l e 2 ) (e m \) (9) 

The general expression ([3]) derived from the score ma- 
trix ([2]) for the average fidelity of the transfer, completed 
at tf, is then (see Suppl.): 



f{tf) = l- dL0G{L0)F{tf,Ul), 



(10) 



=3 



t 2 



dr cos (4>(t)) e" 



dr sin(2^>(r))e 



(11) 



being the transfer-oriented modulation control spectrum. 

Thus, the average fidelity of the transfer has an in- 
volved dependence on the modulation V(t) and the trans- 
fer time tf. The problem at hand is to find the optimal 
transfer that minimizes the average infidelity at time tf, 
1 — F(tf). Obviously, zero infidelity is obtainable for 
infinitely fast (zero-time) transfer, if we allow infinitely 
strong control. Since this is unphysical, we add a con- 
straint on the total energy E of the transfer process 



dt{V{t)f 



E. 



(12) 



As discussed below, this constraint can prevent leakage 
to levels out of the operational qubit subspace [10. The 

constraint defines the minimum possible time for the 
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transfer t mhl = |g. 

We illustrate the general expressions (TlT))) - (rT2"j) for a 
typical non-Markovian Lorentzian bath spectrum, i.e. 
an exponentially decaying correlation function $(t) = 
J. e -|t|/*c j f c being the correlation (memory) time. One 
might expect that for such a simple bath the best strat- 
egy is the fastest possible transfer under the energy con- 
straint, i.e. when the modulation is given by V(0 < t < 
tmin) = 2E/ir. Surprisingly, a slower transfer (tf > t m - ln ) 
with an appropriate modulation <j)(t) (detailed below) can 
improve the average fidelity even for a purely Marko- 
vian bath, with negligible correlation (memory) time 
tc/t m in — > 0, and more so for baths with memory times 
longer than the transfer time, t c > t m { n . 

When the bath is memoryless, i.e. Markovian, this 
improvement is limited, as shown in Fig.[T] to about 12%. 
By comparing the "best" solution to the "fastest" one 
(Fig. ffl(a), (b)), one can see the the "best" solution starts 
off faster and then slows down, being overtaken by the 
"fastest" solution only at t ~ 0.9i m i n . This illustrates 



3 



the source of the Markovian noise reduction: the "best" 
solution starts off faster, so as to transfer more of the 
information while it is still nearly untainted by the bath. 
Obviously, towards the end it must slow down so as to 
comply with the energy constraint, thus resulting in total 
transfer time tf that is longer than the fastest time i m ; n 
for the given energy. 

However, when the memory-time t c of the bath is com- 
parable to or larger than the characteristic transfer time 
tc ^ frmin) a much larger improvement can be achieved 
(see Fig. [I]). Remarkably, the best solution actually per- 
forms a full transfer, <j)(t) = tt/2, well within the mod- 
ulation time, but rather than stop at cj> — tt/2 it then 
"overshoots'" the transfer, so that <j>(t) > tt/2, and then 
returns slowly to tt/2. This can explain the source of 
the noise reduction — when "overshooting" , the infor- 
mation partially returns from the target (storage) qubit 
to the source (noisy) qubit, but with a negative sign. 
Hence, similarly to the "echo" method, the noise now op- 
erates in the reverse direction, correcting itself, i.e., the 
non-Markov bath effect is undone. This requires transfer 
times significantly larger than the minimal transfer time 
£ m in, ranging from 3t m i n to even 10i m i n or more, yet the 
fidelity increases substantially (up to 50% in Fig. [1]). 

Using the Euler-Lagrange variational method one can 
find an analytical solution for the optimal modulation 
phase 4>(t), given a Markovian bath at long times (see 
Suppl. Info.). This yields: 



dx 



' sin 



\24> M {x)) 



2 cos 4 (0 M (x)) 
3 



(13) 



with 0a/(O) = 0. Eq. (|13p determines the shape of 
4>m{x) and its formal "energy" eu = J Wm^- )^^ = 
1.038 . . . (where both x and eu are dimensionless). The 
general Markovian optimal modulation at infinite time 
for any energy E is then <f>(t) = 0jif (— i), with an infi- 



delity of 7^#- = 7 1 -°^ 7 --- ; >y being the dephasing rate of 
the source qubit 1. The fastest modulation with energy E 



has an infidelity of jjrp 



' 8E 7 E 

optimal modulation has about 12% less infidelity than 
the fastest modulation for the same energy. 

Let us now take into account the breakdown of the 
RWA as a noise source. A realistic coupling of the form 
Ox ® cr^\ as in has non-RWA terms of the form 
\g192) {e\e2\ which do not conserve excitation. In all con- 
trol scenarios such terms are either discarded or, at best, 
any transfer of population via non-RWA terms is consid- 
ered alongside all other forms of leakage. However, there 
is a drastic difference in timescale between the non-RWA 
terms in Q and leakage to higher levels: the qubit level 
separation ujq is often orders of magnitude smaller than 
the separation to the next (non-qubit) level. Hence, leak- 
age to higher levels requires timescales that are orders of 
magnitude smaller than those breaking the RWA. 



. This means that the 



0$ 



(&) MarkoviaA bath:' tcAm'in = ' 



/^min — 



(b) Non-Markovian bath: t c /t m - m = 1 - 




tf/t tt 



FIG. 1. Inset - scheme of coupling between "noisy" source 
qubit 1 and quiet target qubit 2: 2-photon transfer off- 
resonantly through \i) gives effective a\ ® a^. coupling with a 
controllable strength V(t) — Kn S^ , where Q(t) is the Rabi fre- 
quency of an external laser field. Main panel: dependence of 
the lowest achievable average infidelity on the transfer time t / 
normalized to the fastest transfer time i m i n at a given trans- 
fer energy (Eg. I12|) . This function is plotted for various bath 
memory times: (black, solid) t c — (Markovian); (red, dash) 
t c = tmin', (green, dash-dot) t c — 10i m m. Even for Markovian 
baths (t c <C t m in) the best solution is not the fastest one. 
For non-Markovian baths (t c > t m j n ) two plateaux (regions 
of insensitivity to tf) can be seen. The first plateau is in- 
dependent of the memory time, and matches the Markovian 
plateau. The second plateau is lower the longer the memory 
of the bath. In (a) and (b) the transfer phase <j>(t) is plotted 
versus t/t m in. The fastest modulation (black, dotted) with 
the Markovian optimal modulation (red, solid) and the non- 
Markovian optimal modulation (green, dashed), in Marko- 
vian (a) and non-Markovian (b) baths. In the Markovian 
bath the optimal modulation transfer starts off faster than 
the "fastest" transfer (when the information is still "fresh"), 
and slows down subsequently. For the non-Markovian bath, 
optimal modulation achieves full transfer (<j>(t) = tt/2) well 
within the modulation, but then "overshoots" (4>(t) > tt/2), 
and eventually returns to <j>(t) = tt/2. 



The RWA in (T5)) breaks down if the transfer time is 
similar to or less than the inverse energy separation of 
the qubit, t m j n < ^q 1 - This is a common case for qubits 
whose resonance frequency is in the microwave (GHz) or 
RF (MHz) range. In such cases, the optimization process 
must take into account both the dephasing due to the 
bath as in (fl~0)) - (fl"2"l) and the error due to the non-RWA 
terms when minimizing the infidelity. 

If one of these errors is vastly larger than the other it 
is the only one to consider. The problem changes only 
if the bath-induced and the non-RWA errors are similar. 
In this case we find that dephasing of the doubly excited 
state caused by erf (g> B in ([6]) is a fourth order effect and 
hence can be ignored in the present second order treat- 
ment. The result of this approximation is that the system 
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FIG. 2. Fastest (black, dots) vs. best modulation for RWA 
(red, dash) and non-RWA (green) transfer with a final time 
of tf — 5t m in (top) and tf = 10t m in (bottom) in a non- 
Markovian bath. The level separation ojo satisfies wotmtn = 7T, 
giving ~ 2.5% population leakage for fastest modulation and 
~ 10% loss of fidelity from decoherence. The sinusoidal "wig- 
gling" in the non-RWA solution is larger when the final time 
is shorter (top). 



can be split into two completely separate subsystems O 
and £, the former suffering only from dephasing and the 
latter only from unwanted population of the doubly ex- 
cited level \e\e2). The Hamiltonians of these systems are 



o 



B, 



(14) 



The goal of our optimization is to find a control V(t), 
shared by both subsystems, which maximizes the fidelity 
of the transfer in subsystem O (as per (|10 p -(|12 J> ) while at 
the same time minimizing the doubly-excitation in sub- 
system £. 

The optimal modulation for Markovian and non- 
Markovian baths, for different final times tf, is given in 
Fig. [2j The result shows that the optimal modulation 
resembles the solution for dephasing in Fig. [1] but with 
added "wiggles" . The "wiggling" takes up the entire time 
allowed for the transfer, but, given more time for the to- 
tal transfer, the amplitude of the "wiggle" diminish. This 
can be understood as follows: first you should complete 
the transfer assuming the RWA so as to minimize the in- 
formation lost to the bath. Once the transfer is complete, 
and decoherence is minimized, we can use whatever en- 
ergy is left to return the "leaked" excitation from |eiea) 
back to Iffiffa)- This can be done by a weak sinusoidal 
modulation of frequency 2u>o, inducing a Rabi coupling 
between the doubly excited and zero-excited levels of £ . 
In practice this is not a two-stage modulation, as the 
weak oscillation is superimposed on top of the transfer 
modulation. 



The energy needed to "undo" the non-RWA effect is 
inversely proportional to the allowed time — E » 
(where ip ee is the amplitude of the doubly excited state 
|eie2)). Hence, given enough time, the correction of the 
non-RWA effects requires negligible energy, yielding the 
same results as in the RWA case. If, however, time is 
limited — a larger proportion of the energy of the trans- 
fer must be reserved for the correction of the non-RWA 
effect, resulting in smaller reduction of the bath-induced 
noise. 

To conclude, our analysis of state-transfer optimization 
within hybrid open systems, from a "noisy" qubit to its 
"quiet" counterpart, has revealed an intriguing interplay 
between our ability to avoid both bath-induced errors 
that profoundly depend on the bath-memory time and 
the limitations imposed by leakage out of the operational 
subspace. Counterintuitively, under no circumstances is 
the fastest transfer optimal (for a given transfer energy). 
Generalizations to higher-dimensional cases are expected 
to follow analogous trends. 
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